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ABSTRACT 

An attempt to measure the Hubble constant with gravitational lens time delays is often 
limited by the strong degeneracy between radial mass profiles of lens galaxies and the 
Hubble constant. We show that strong gravitational leasing of type la supernovae 
breaks this degeneracy; the standard candle nature of type la supernova luminosity 
function allows one to measure the magnification factor directly, and this information 
is essential to constrain radial mass profiles and the Hubble constant separately. Our 
numerical simulation demonstrates that the Hubble constant can be determined with 
^ 5% accuracy from only several lens events if magnification factors are used as 
constraints. Therefore, distant supernova survey is a promising way to measure the 
global Hubble constant independently with the local estimates. 

Key -words: cosmology: theory — distance scale — galaxies: structure — gravita¬ 
tional leasing — supernovae: general 


1 INTRODUCTION 

Gravitational lensing is known to be a powerful tool to de¬ 
termine the Hub ble constant H p directly, without using the 
distance ladder ( R,efsda| 1964). Derived values of Hq show 
good agreement among known several lens syst ems, once the 
radial mass pro fi les of lens galaxie s are fixed (Koopmans & 
[Fassnachi 199!;; Kochanek 2002a). If the lens galaxies are 
assumed to have the singular isothermal mass distribution, 
analysis of five gravitational lens syst ems indicates tha t the 
value of Ho is Hq ~ 50kms“^Mpc“^ (Kochanek 2002a) and 


hence too low to be cons istent with the local m easurement 
Ho ~ 70kms“^Mpc“^ ( Freedman et al. 2001 ). This dis¬ 
crepancy niay be ascribe to inhomogeneity in the universe 
(Pomita 2000a and therefore it is important to study the 
global Hubble constant independently with the local mea¬ 
surement which relies on the distance ladder. 

The main limitation of this technique is that there is 
a strong degeneracy between radial mass profiles and Hp 


( WambtjRanHH & Pac^ynski ]J194t_ i^tofi & Kodmnek| 19971 

Koo£m M^fc JjAssnadit - I-99S - Wittj^laOj fcJ£eeto^ 


2002t [Wucknitz ^002| ; [Kochaii^ 2002a; Zhao & 


Tada & Futamase 200C p Wiliiams_&_Sal^ 20001 Chiba & 


Takahashi 

Qinf 200^ Therefore, unless we specify the radial mass dis- 
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tribution in the lens object, we hardly constrain the value 
of Ho- The strong dependence of differential time delays on 


the radial mass 
that statistics o 
density profiles ( 
broken if the Eii 

distribution 
time delay 
Oguri et al. 
istein ring 

, on the other hand, indicates 
s provide a powerful probe of 
2002). This degeneracy can be 
mages of host galaxies are ob- 

served ( 

Kochanek, Keeton, & McLeod 

2001 

), but the obser- 


vation of host galaxies is often very difficult because of the 
large brightness contrast between quasars and host galax¬ 
ies. Other way to br eak this degenerac y comes from the 


central core image s (Rusin fc Ma 


2001 


Evans & Hunter 


2002 ; Keeton 200f). The lack of central core images, how¬ 
ever, places only the lower limit of the mass concentration, 
and this corresponds to the lower limit of Hq. Information of 
stellar kinematics and the mass-to-light ratio also allows one 
to constrain the radial mass profile of the lens galaxy and to 
break the degeneracy in the Hubble constant, although the 
measurement of velocity dispersions for dist ant galaxies is 
very d i fficul t and involves large uncertainties. Treu & Koop- 
mant (2002|), for instance, concluded that the mass density 


profile of the lens PG1115-I-080 is steeper than the singular 
isothermal mass distribution, and obtained a value of the 
Hubble constant Ho ~ 60kms“^Mpc“^ which is marginally 
consistent with the local measurement. 


In this Letter, we consider strong lensing of distant su¬ 
pernovae (SNe) whi ch is expected to be observed by future 
observational plans (Kolatt & Bartelmann 1998; Wang 200C; 


© 2002 RAS 

























































































L2 M. Oguri and Y. Kawano 


Porciani & MadarJ tOOOl; Sullivan et al. |200C| 

IHoIz 

2001 

; Goo- 


the satellite (SuperNova/Acceleration Probe) will 

catch at least a few lensed SNe per year. The well-known 
advantages of gravitational leasing for distant supernovae 
include: ( 1 ) differential time delays between images can be 
easily measured with high accuracy (of order of one hour), in 
marked contrast with usual quasar leasing for which time de- 


lay measurement is quite difficult (e.g., Kundic et al. 1997); 


(2) since SNe are transient phenomena, it is easier to ob¬ 
serve lensed host galaxies of SNe after SNe are faded away. 
More importantly, we find that information of magnifica¬ 
tion factors, which can be directly measured in the case 
of type la SN leasing, breaks the degeneracy between ra¬ 
dial mass profiles and Hq. First we provide the reason for 
this degeneracy breaking analytically, and next we numeri¬ 
cally demonstrate the importance of magnification factors. 
In what follows, we adopt a lambda-dominated universe with 
(flo, -^o) = (0.3, 0.7), where is the density parameter and 
Ao is the cosmological constant. The Hubble constant in 
units of 100 kms“^Mpc“^ is denoted by h. 


2 ANALYTIC CONSIDERATION 


In this section, we estimate how accurately we can constrain 
radial mass profiles and the Hubble constant separat ely, on 
the b asis o f simple analytic consideration. Following Wuck- 
nitz (M 02 ), we consider the degeneracy between radial mass 
profiles and the Hubbl e constant in terms of well-know n 
mass-sheet degeneracy ( Falco, Gorenstein, fc Shapiro 1985 ). 
First start from the lens equation: 


y = Xi- 


( 1 ) 


where y and Xi are the positions of the s ource and images 
and ^(xj) de notes the lens potential (see 


Schneider, Ehlers 


&: Falco 1992). If we multiply the lens equation with (1 — k) 
then the lens equation can be rewritten as 


(1 - K)y = Xi-V 


(1 - K)tp{Xi) + K-^ 


( 2 ) 


Therefore, the image position Xi is never changed if we 
transform unobservable quantities as y —> (1 — K)y and 
ip{xi) ^ (1 — K)'ip{xi) + K{xif2). This is the mass-sheet 
degeneracy. Since the time delay is calculated from 

hAti oc -(3) 

this transform also changes the estimation of the Hubble 
constant, /i —> (1 — k:)/i if Atij is fixed to the observed value. 
Therefore, the Hubble constant h cannot be uniquely deter¬ 
mined from information of {xi} and {Atij} = {Ati — Atj}. 
In the usual quasar leasing, additional information is sup¬ 
plied by the flux ratio. The flux ratio is simply derived from 
the ratio of magnification factors yt: 


m = 


dxi 


(4) 


From this expression, it is found that y,i is transformed 
as Hi —> {1 — k)~^ fii. Hence the flux ratio rij = yt/yj is 


^ http://snap.lbl.gov/ 


also never changed. This means that information of nj can¬ 
not break this degeneracy. The degeneracy between radial 
mass profiles, tp (x , and the Hubble constant h can be 
interpreted in this context b ecause of the simple relation 
1 — K = 2 — P (Wucknitz 2002). This yields a general scaling 
law, 


/i oc 2 — /3, 


(5) 


without changing observable values such as Xi, Atij, and 


In the case of SN la leasing, the situation changes dras¬ 
tically. “Standard candle ” nature of SNe la (Phillips 199f; 
Riess, Press, & Kirshnei 1996) allows one to observe mag¬ 


nification factors yi directly. The magnification factor has 
the strong dependence on the r adial mass profile (see also 
Wambsganss & Paczynski 1994; Kawano et al., in prepara¬ 


tion), 

Pi oc (2 - P)~ 


( 6 ) 


thus we can break the degeneracy b y the measurement of 
magni fication factors. We note that Kolatt & Bartelmann 
( [l998| ) also proposed this method to break the mass-sheet 
degeneracy in galaxy clusters and reconstruct galaxy cluster 
mass. Although yi should have significant error which arises 
from the intrinsic dispersion in SNe la pe ak luminosities as 
well as substructure in the lens galaxy (Mao fc Schneider 


199S ), its effect on h estimation is not so severe because 
and 


equations 
a half of that of yi. 


implies that the resulting error of h is 


3 SIMULATED RESULT 


To illustrate how accurately we can determine /3 and h, we 
show the results of our simulation. First, we assume the lens 
galaxy is well characterized by the following lens potential 
(Kawano et al., in preparation): 

'ip{x) = r^(ao-I -02 cos 20-I -03 cos 30-I -62 sin 20-I -63 sin 30), (7) 


where x = (r, 0 ) is the position in the polar coordinate. 
We note that best-fit values of h and y in the analysis of 
quasar PG1115-I-080 are quite ins ensitive to the cho ice of the 
angular part of the lens potential (Kochanek 2002b| ; Kawano 
et al., in preparation), thus it is sufficient to analyze this lens 
potential only. We also take account of the effect of external 
shear as 


^shear(a;) = ^7’’^ cos (20 - 20 ^). 


( 8 ) 


We randomly put a source, generate quadruple images, and 
calculate differential time delays and magnification factors, 
assuming following parameters: /3 = 1.0, h — 0.5, ao = 0.5, 
02 = b 2 = 0 . 01 , as = bs = 0 . 001 , 71 = 7 cos 20 -,. = 0 . 1 , 
and 72 = 7 sin 20 -, = — 0 . 01 , where at and bi are values 
when X is in units of arcsec. Source and lens are placed at 
as = 1.5 and zl = 0.5, respectively. Generated images have 
separations on the order of . For the observable quantities, 
such as image positions and time delays, the Gaussian noise 
is added. We assume following dispersions: 


a, = 0.01", 

(JAt = 0.05 [day], 

O'logr^j — 0.09, 


(9) 

( 10 ) 

( 11 ) 
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Figure 1. Constraints on the radial mass profile /3 (eq. [Q) and 
the Hubble constant h. The contours of Ax^ = ~ Xmin 

l3-h plane are calculated from one quadruple lens event. Crosses 
indicate the assumed value in generating observable quantities; 
{f3,h) = (1.0, 0.5). Upper panel, image positions Xi, differential 
time delays and magnification ratios rij are used to calcu¬ 

late x^ • Lower panel instead of rij , magnification factors pi are 
used to calculate x^- Dotted lines are same as solid lines, but in 
this case additional non-Gaussian errors due to microlensing are 
also included. 


O'log fli 


= 0 . 12 . 


( 12 ) 


The precisions in positions and time delays are consistent 
with 0.1 pixel of the instrument and estimated accu racy for 
SNe la ligh tcurves in SNAP survey, respectively (Goobai 


et al. i002). The dispersion of magnification ratio, which 
roughly corresponds to ~ 20% fractional error, is a fidu¬ 
cial er ror often assumed in x minimization (e.g., Kochanek 


2002a). We assume that the dispersion of the magnification 
factor is somewhat larger than this, roughly corresponds to 
~ 30% fracti onal error, because not o nly substructure in the 
lens galaxy (Mao & Schneider 1998) but also the intrinsic 
dispersion of SNe la peak magnitudes contribute to criog^i;. 
Other possible source of the dispersion is dust extinction in 
the lens galaxy. However, the effect of dust extinction can be 
corrected from the observed redde ning because of kn owledge 
of an SN la’s intrinsic color (e.g., Riess et al. 1996). 


After the virtual “observational data” is generated, we 
perform minimization using the same lens model. At that 
time, we fix values of /3 and h, and optimize the other pa¬ 
rameters such as fli, bi, 7, and the source position. Faint 
core images which may appear when /3 > 1 are always ne- 
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Figure 2. Same as Figured, but from five quadruple lens events. 


glected. In calculating we consider following two cases: 
(1) Only the magnification ratio r is measured. This case 
corresponds to traditional quasar lensing. (2) The magnifi¬ 
cation factor is directly measured. This is the case of SN 
la lensing we are interested in. For each case, we calculate 
the contour of — Xmin in fhe P-h plane. Figure ^ 

plots constraints on /3 and h from one quadruple lens event. 
This figure clearly shows that in the case of SN la lens¬ 
ing P and h are well constrained separately. It is surprising 
that the Hubble constant is determined with ~ 10% accu¬ 
racy (68% confidence) from only one lens system. On the 
other hand, when magnification factors are not used, P and 
h are poorly determined; they show the strong degeneracy 
h oc 2 — p. We note that in practice constraints from quasar 
lensing may be worse than our result using magnification ra¬ 
tios, because the error of time delays is usually much larger 
than our assumption (eq. |^). Figure ^ shows constraints 
from five quadruple lens events. In generating observables 
for each event, the position of the source is changed while 
the lens model is always fixed. In this figure, the Hubble con¬ 
stant h is determined with ~ 5% accuracy (68% confidence) 
when magnification factors are used, while the accuracy is 
still ~ 20% (68% confidence) when magnification ratios are 
used. We also examine the case that lens galaxies have dif¬ 
ferent values of P, and the result is shown in Figure ^ In this 
plot, we assume that five lens systems have different radial 
mass profiles; P = 0.8, 0.9, 1.0, 1.1, and 1.2, respectively. 
This figure clearly indicates that the magnification factor 
is quite useful to constrain the Hubble constant even if the 
scatter of P is taken into account. The contour is slightly 
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Figure 3. Same as Figure but from five quadruple lens events 
with different radial mass profiles; (3 = 0.8, 0.9, 1.0, 1.1, and 1.2. 


elongated along the /3-direction by the scatter of /3, but the 
accuracy of the Hubble constant determination is almost 
same as that in Figure ^ On the other hand, the strong 
degeneracy still remains when only the magnification ratio 
is used as constraints. Therefore, we conclude that magnifi¬ 
cation factors which are observed in SN la leasing provide 
indeed important information to break the /3-/i degeneracy. 

The most important effect we neglect here is the effect 
of microlensing by stars in the lens galaxy. The magnifica¬ 
tion probability distribution of microlensing has a long tail 
at high magnification regions and thus induce non-Gaussian 
errors. This means that microlensing may change our re¬ 
sult systematically. The robust estimate of the magnification 
probability distribution of microlensing is, however, difficult 
because the amount of microlensing fluctuation depends on 
the fr action of the stellar componen t which is still unclear 
(e.g., Bchechter & Wambsganss 2002). Therefore, in this pa¬ 
per we show the effect of microlensing by simply adding 
non-Gaussian errors. We assume the following log-normal 
probability distribution: 


P(log/rML) = 


exp {-[ln(log/rML + m)]^/2s^} 
V2ns^{\og fiMh + m) 


(13) 


where /tml indicates the additional magnification due to mi¬ 
crolensing. The parameters m = 1 and s = 0.2 are chosen so 
as to reproduce the typical fraction which suffers from mi¬ 
crolensing variability; in our model the fraction magnified 
by log/TML > 0.4 (> 0.2) is ~ 5% (~ 20%) respectively, and 
this fraction is roughly consistent with the previous estima¬ 


tions (^oopmans fc Wambsganss 2001; Wyithe & Turner 


2002). The probability distribution of microlensing should 
depend on convergence and shear at images and therefore 
each image may have different probability distributions. We 
neglect this effect because we now consider quadruple lens 
systems in which images have approximately the same sep¬ 
arations from the lens center and are likely to have similar 
convergence and shear. Gonstraints on /3 and h including 
the effect of microlensing are shown in all figures (dotted 
lines). These figures indicate that the effect of microlensing 
is fairly small. The magnification by microlensing systemat¬ 
ically moves contours mainly to the lower h, but the devi¬ 
ations are sufficiently small and our assumed model (/3 = 1 
and h = 0.5) lies still within contours. Our conclusion is 
therefore that the effect of microlensing is not so severe. 


4 SUMMARY AND DISCUSSION 


We have shown that the strong degeneracy between radial 
mass profiles j3 and the Hubble constant h can be broken 
if we use magnification factors as constraints. This means 
that SN la leasing has the great advantage over traditional 
quasar leasing. We have found that in the case of SN la 
leasing the Hubble constant is constrained quite accurately, 
with ~ 5% accuracy from only several lens events. We have 
found also that both the scatter of radial mass profiles and 
microlensing do not affect our result so much. In contrast 
to this, quasar leasing can poorly constrain the Hubble con¬ 
stant, even if the same accuracy of time delay measurements 
is assumed. Of course, this method can be applicable to any 
distant astronomical objects which have quite narrow lu¬ 
minosity function. The limitation of our method is that it 
probes only the local slope of the mass profile. Therefore, 
the value of the Hubble constant derived from our method 
may be different from the true value if the radial mass profile 
is significantly different from a power-law. 

Although SN leasing has not ever been observed, lensed 
SN will be found in the future observations such as SNAP. 
Then how many lensed SNe are expected to be observed in 
the future observational plans? SNAP survey, for instance, 
can catch 2000 ~ 3000 Type la SNe at z < 1.7 per year. 
Since the leasing probability at 2 ~ 1.5 i s ~ 10~^, we ob¬ 


serve at lea st a few lensed SNe per year (Holz 2001; Oguri 


et al. 


2003). Time delays between images are always ob¬ 


served in the case of SN leasing because SNe are transient 
phenomena. Hence all lensed SNe can be used to constraint 
the Hubble constant. Therefore, a few year’s observation by 
SNAP is enough to derive the accurate value of the Hubble 
constant. Note that our method can constrain the Hubble 
constant accurately from even one lens system. On the other 
hand, large-scale surveys, such as two degree Field system 
(2dF) and Sloan Digital Sky Survey (SDSS), will also find 
more than one hundred of quasar lenses. The number of 
lens systems for which time delays are measured, however, 
should be much smaller than this because of the difficulty in 
measuring time delays. Moreover, additional observations, 
such as host galaxies of quasars or velocity dispersions of 
the lens galaxies, are needed to break the /3-/i degeneracy. 
Therefore we conclude that our method is practically appli¬ 
cable and indeed have advantages compared with methods 
using quasar lenses. 
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In our simulation, we considered only quadruple lens 
systems which have larger numbers of constraints. In prac¬ 
tice, SNAP-like survey will catch double lenses much more 
than quadruple lenses because the magnihcation bias is ne¬ 
glected at least for \aw-z SNe {z < 1.7 in the case of SNAP 
survey). Time delay bias (Oguri et al. 2003) favors quadruple 
lenses, but the time delay bias is not so significant at typi¬ 
cal lens separation 0 ~ 1”. We believe, however, our results 
are also applicable to double lenses, because magnihcation 
factors are sensitive to radial mass prohles even when the 
source is far from the center of the lens and thus likely to 
be the double lens (see Oguri et al. 2002). 


Finally, we comment on the measurement of the mag¬ 
nihcation factor /i. Consider the situation that only lensed 
SNe are observed and the absolute magnitude of (unlensed) 
SNe, w hich is estimated using the local Hubble constant, is 


known. In this case, the hux of unlensed SNe at the same 


redshift, expected from the absolute magnitude, scales as 
/unlensed oc if the assumptiou that the local Hubble con¬ 
stant is the same as the global Hubble constant (denoted by 
h) is relaxed. Since we observe /lensed directly, the magnihca¬ 
tion factor should scale as /r = /lensed/Zunlensed oc h~^. This 
scaling is exactly same as that derived in this means 
that the f3-h degeneracy is never broken in this situation. 
But if the redshift of lensed SN is not so larger, huge num¬ 
bers of unlensed SNe which have the similar redshift should 
be also observed. In this case, /r can be simply estimated 
from the magnitude diherence between lensed and unlensed 
SNe independently with h. The only assumption needed to 
estimate fj, is therefore that SN la is an excellent standard 
candle. We emphasis that several important uncertainties 
in the SN survey, such as cosmological parameters, possible 
evolution of intrinsic luminosities, dust extinction outside 
the lens galaxy, and the Cepheid calibration, do not affect 
the measurement of p. 
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